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Abstract 

We present a Maplell+GRTensorll based symbolic calculator for 
instanton metrics using Newman-Penrose formalism. Gravitational 
instantons are exact solutions of Einstein's vacuum field equations 
with Euclidean signature. The Newman-Penrose formalism, which 
supplies a toolbox for studying the exact solutions of Einstein's field 
equations, was adopted to the instanton case and our code translates 
it for the computational use. 

Keywords : Gravitational instantons; Newman-Penrose formalism; Symbolic 
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1 Introduction 



The interest in symbolic computational study of general relativity is growing 
rapidly as the capacity of the the computer systems increase. The computer 
is no longer an apparatus for the numerical relativist, but as the symbolic ma- 
nipulators get easier to use, the more researchers get into the subject by using 
these systems. Comprehensive historical reviews can be found in references 
[H [21 [3]. Another important point to consider is the specialized packages 
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which run under these platforms. GRTensorll pE], being a widespread pack- 
age in general relativity led significant progress in the area El El El El EE] • 
The Newman- Penrose (NP) formalism supplies a toolbox for investigating the 
exact solutions of the Einstein's field equations [TU [12] . Goldblatt has devel- 
oped NP formalism for gravitational instantons based on the SU(2) x SU(2) 
spin structure of positive definite metrics [T51 HI) . In the gravitational in- 
stanton case, the gravitational field decomposes into its self-dual and anti- 
self-dual parts and this decomposition is natural in the spinor approach 
which necessitates two independent spin frames for the spinor structure of 
4-dimensional Rieman manifolds with Euclidean signature [15] . 
Gravitational instantons are exact solutions of Einstein's vacuum field equa- 
tions with Euclidean signature. They are analogous to the Yang-Mills instan- 
tons, the finite action solutions of the classical Yang-Mills equations. They 
admit hyper-Kahler structure [T5| [16]. For a detailed overview one can see 
Eguchi et al's review and the articles cited in this paper [PT] . 
Aliev and Nutku applied differential forms to the NP formalism for gravita- 
tional instantons which made the formalism more suitable for the symbolic 
computation [T5] . 

2 The program 

Our program, NPInstanton is designed for calculating physical and mathe- 
matical quantities for instanton metrics using Newman- Penrose formalism. 
It is coded under Maple 11 and GRTensorll package. The quantities which 
can be calculated are: 

- massless scalar equation, 

- massless Dirac equation, 

- source-free Maxwell equations, 

- covariant and contravariant Dirac 7 matrices, 

- coframe /(= l^dx^ 1 ) and m(= m fl dx fl ), 

- Ricci rotation coefficients, 

- Weyl scalar s, 

- trace-free Ricci scalars, 

- spinor equivalent of the connection 1-forms, 

- basis 2-forms, 

- curvature 2-forms, 
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- integrands of the Euler number and the Hirzebruch signature curvature 
part integrals, 

- Petrov class of the space. 

As one can see, some of the objects could be calculated by standard means 
and without using a signature-dependent package. But for the sake of com- 
pleteness, we added these features to the program. By these, the program 
becomes a complete symbolic calculator for an instanton metric. The NP cal- 
culator of GRTensorll is not designed for the Euclidean signature and could 
give unexpected results. Therefore, a complete NP based calculator, com- 
bining the power of GRTensorll and NP formalism for these special metrics 
is useful. 

The program uses, 

1. GRTensorll Package: Several objects (Ricci scalar, covariant Weyl tensor, 
etc.) are calculated by this package in our program. No specific knowledge 
of this package is needed for using NPInstanton. Our program creates the 
metric file needed for the calculation itself using the NP legs given in the input 
file and writes it to the metric directory of GRTensorll as "npinstanton.mpl". 
This metric file then can be used in GRTensorll independently. 

2. DifferentialGeometry Package of Maple 11: The definition of the wedge 
product is supplied by this package in our program. The wedge product 
definition is also supplied by "difforms" package in Maple but it has conflicts 
with linear algebraic quantities. This package is excluded in a limited version 
of the program for the users of older Maple versions. 

3. linalg package of Maple: This rather old internal package is used for an 
eigenvalue calculation in Petrov classification section. 

The program is set for a computer having a 512 Mb of RAM (Average value 
for today's personal computers). If the system has less memory, the user 
must change the line 

kernelopts (gcf req=10~7) : 

of the npinstanton.mws file to a lower value (KT6 is the standard value 
of Maple). For a computer having larger than 1 Gb of memory, the user 
may change the gcfreq value as KT8. "gc" is the abbreviation of "garbage 
collection" and it is the Maple's internal routine which cleans the memory 
after an amount of memory is allocated. For a computer having a large 
amount of memory, one can increase the frequency of this process. The 
larger the gc frequency value results in more memory to be wasted but for 



3 



a system having a large amount of memory it increases the performance for 
some calculations. 

The programming style is procedure-based. Each command is a Maple pro- 
cedure which can call other calculation procedures and inherit their outputs. 
For example, when the dirac() command is given for the calculation of mass- 
less Dirac equation, the program calls dirac() procedure. The calculation 
of Dirac equations needs the information on 7 matrices. Therefore, dirac () 
procedure calls gammamatricespless() procedure for this calculation and 
its output is sent to dirac () procedure. The suffix "-pless" (printless) means 
that this procedure is a "secondary" one and it does not print its output to 
the screen. When dirac () procedure finishes execution, it removes the in- 
formation that was inherited from gammamatricespless() procedure and 
sends the user only its permitted output as dirac vector. 
The "secondary" procedures do the massive works of the main procedures 
and the main procedure generally contains the key directives. Therefore, this 
usage makes the program easier to track. In addition to this, if a user wants 
to add another command to the program, a procedure of that command can 
be coded using these "secondary" procedures, or the user can add his/her 
own "secondary" procedures, as well. 

Throughout the program, only the output variables of the procedures are the 
global ones and only these variables can be reached by the user at the end of 
the execution of the command. Therefore, the user should be aware of the 
output names of the procedures that are being used. This global variables 
are shown at the end of the execution of each command. For example, the 
output variable of scalaroperator() command will be shown as scalarop 
and it can be called in the session when needed. 

The outputs of the auxiliary calculations are cleared from the memory by 
setting them as local variables, as these auxiliary calculations may occupy 
too much memory. For example, 7 matrices information is not necessary if 
the user needs to know only the Dirac equation. Therefore, dirac () proce- 
dure clears the information about 7 matrices after finishing the execution. If 
the user needs to know the 7 matrices, gammamatrices() procedure is run 
by the user. As another example, in connlform() procedure, only output 
values are connection 1-forms and the outputs from coframepless() and ric- 
cirotcoeffplessQ procedures, as well as internal calculations are dismissed 
to gain memory for further calculations. 

A few commands for simplification are added to the code which can be eval- 
uated easily by an average personal computer (i.e. having a 512 Mb of 
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memory), but it is always more convenient for the user to choose the right 
simplification technique for the problem after the calculation. The output 
must be regarded as a "raw material" for a simplification routine that is to 
be chosen by the user. The user having a computer with insufficient memory 
can extract the simplification routines from the program, simply by modify- 
ing it in an editor but it is not recommended as it may lead to miscalculation 
of some properties such as Petrov type or (anti-) self-duality. 
One can reach the program files using the web site given in the "Final Re- 
marks" section. Linux and Windows versions are available and for those who 
have older Maple versions, a limited version of NPInstanton which does not 
contain the basis2form(), curv2form() and topologicalnumbers() parts 
(parts that use DifferentialGeometry Package) is also supplied. The details 
on how to run the program is given in the README file which comes with 
the program files. 



2.1 The input file 

The input file contains the necessary definitions for the space and wave equa- 
tions. General definitions (constants, etc.) can also be given in this file. The 
definitions should be given using Maple's syntax. The input file of our pro- 
gram plays the role of the metric file of GRTensorll. 
The user first sets the coordinate names in the first section as, 

f irstcoordinate : = x: 
secondcoordinate : = theta: 

etc. as given in example input file. Then the components of the covariant NP 
legs are entered. l_covar[l] being the first component of the covariant Z M leg 
and l_bar_covar[l] being the complex conjugate of l_covar[l], etc.. Maple 
does not do the complex simplifications because they need assumptions that 
may cause wrong calculations. Therefore, the most appropriate way to define 
the legs is to give both by hand. The choice of NP legs is not unique for a 
metric. For the Eguchi-Hanson NP legs we can have[17j: 
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Then the input forms are 

l_covar[l] :=l/(sqrt(2)*A) : 
l_bar_covar[l] : =1/ (sqrt (2) *A) : 

etc. as given in example input file. The A value being \Jl — ^ was given in 
the general definitions section. 

Spinor components can be given if the user will be calculating the Dirac 
equation. They can be chosen as 





_ e K mJ < 
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(3) 
(4) 
(5) 
(6) 

and the input form is 

spinor component 1 : =exp(I* (m+(l/2) ) *xi) *psil [r , theta,phi] : 

etc. as given in example input file. 

The scalar function can be given for the calculation of the scalar equation. 
This definition can be skipped if the user does not need to calculate this 
object. It can be chosen as 

y? = e im «$(r,M) (7) 
for the Eguchi-Hanson space, <p being the scalar function and it is set as 
scalarf unction: =exp(I*m*xi) *Phi (r ,theta,phi) : 
in the input file. 
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2.2 Command definitions 



The list of commands can be given as, 

scalaroperatorO 
diracO 
maxwell () 
gammamat r i c e s ( ) 
cof rame() 
spinrotcoef f () 
weylscalar () 
tf ricciscalar () 
connlf orm() 
basis2f orm() 
curv2f orm() 
topologicalnumbers () 

and the definitions of these commands are the following, 
• scalaroperator() : 

This command calculates the massless scalar equation, finding the scalar 
operator. The scalar operator is given as 



Here, g is the determinant of the metric. 

The procedure takes the scalar function (name: scalarfunction) from the 
input file. 

The output to be used thereafter: 
> scalarop; 

For the massive case, one can equate this object to M 2 (p 2 , M being the 
mass of the scalar particle. As an additional property, this procedure is not 
dependent on the metric signature and it does not use the NP objects so 
it can be used for any space in four dimensions by extracting it from the 
program. 




(8) 



• dirac(): 
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This command calculates the massless Dirac equations, 



7' 



M V M ^ = 



(9) 



where 



(10) 



and are the spin connections. The procedure takes the spinor vector com- 
ponents from the input file under these names: spinorcomponentl (tfii), 
spinorcomponent2 (1P2), spinorcomponent3 (^3) and spinorcompo- 
nent4 (^4). The output to be used thereafter is the components of the 
"dirac" vector as (i=l, 2, 3, 4): 

> dirac [i] ; 

For the massive case, one can equate these objects to ^-ifj vector, M being 
the spinor mass and i = y/—l. 

• maxwellQ: 

This command calculates the source-free Maxwell equations using the equa- 
tions 102-109 of [15] . Fi/s in the output are the usual Maxwell field matrix 
components. The output is set to be "maxwell" vector whose components 
are the source-free Maxwell equations as (i=l, 2, 3, 4): 

> maxwell [i] ; 

• gammamatrices(): 

This command finds the covariant and contravariant Dirac 7 matrices and 
tests them with anticommutation relations 



The output is "gamma_contr" (contravariant 7 matrices: 7^) and " gamma_covar" 

(covariant 7 matrices: 7^) vectors as (i=l, 2, 3, 4): 

> gamma_contr [i] ; 

> gamma_covar [i] ; 



{7^,7*} = 2g- 



(11) 
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The calculation of the Dirac 7 matrices could be put into the dirac() proce- 
dure but for further calculations involving higher dimensions one may need 
only these matrices. To find the 7 matrix for the extra dimension, one can 
use these results and look for the 7 matrix of the extra dimension using the 
anticommutation relations as in reference [19J. Therefore the procedure is a 
separate one. 

• coframeQ: 

This command calculates the coframe I =L (= l^dx^ 1 ) and m =M (= m^dx^). 
The output variables are the following, "_bar" denoting the complex conju- 
gate: 

> L; 

> M; 

> L_bar; 

> M_bar; 

• riccirotcoeffQ: 

This procedure calculates the Ricci rotation coefficients using equation 20 of 
[15] . The output can be used later by calling: 

> npkappa; > nptau; > npsigma; 

> nprho; > nppi; > npnu; 

> npmu; > nplambda; > npgamma; 

> npepsilon; > npalpha; > npbeta; 

in the meaning of the Ricci rotation coefficients k, t, a, p, n, u, fi, A, 7, e, a, 
(5 respectively. 

• weylscalar(): 

This command calculates the Weyl scalars using equation 68 of [15]. The 
output variables are: 

> weylscalarO; > weylscalarl; > weylscalar2, > weylscalar3; 

> weylscalar4; > weylscalartildeO; >weylscalartildel ; 

> weylscalartilde2; > weylscalartilde3 ; > weylscalartilde4; 
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for $ , ^1, $2, ^3, ^4, ^o, ^i, *2, ^4 respectively, "tilde" means the 
variable has a tilde in the meaning of a different spin frame. 
Another property of this procedure is that it finds out the Petrov class of 
the space according to equations 120-122 of [15]. Petrov classification of 
Euclidean spaces were first studied by Hacyan [21] and then by Karlhede 

• tfricciscalar(): 

This command calculates the trace-free Ricci scalars using equation 69 of 
[15] and sets them to variables for a later use: 

> tf ricciscalarOO ; > tf ricciscalarOl ; > tf ricciscalar02 ; 

> tf ricciscalarlO ; > tf ricciscalarll ; > tf ricciscalarl2 ; 

> tf ricciscalar20 ; > tf ricciscalar21 ; > tf ricciscalar22 ; 

> scalofcurv; 

for $ 00 , $oi, $02, $io, $n, $12, $20, $21, $22, scalar of curvature A respec- 
tively. 

• connlform(): 

This procedure calculates the spinor equivalent of the connection 1-forms 
given by equations 36-37 of [15J. The output can be reached by calling 

> GAMMA00; > GAMMA01; > GAMMA10; > GAMMA11; 

> GAMMAtildeOprOpr; > GAMMAtildeOprlpr ; 

> GAMMAtildelprOpr; > GAMMAtildelprlpr ; 

for Tf], Tq, r^, Y\, T%, Tq,, Ty, Ty respectively where "tilde" means the 
variable has a tilde and "pr" means "prime" in the meaning of a different 
spin frame. 

Another property of the procedure is that, it finds out the gauge by checking 
the necessary and sufficient conditions for (anti-)self-duality namely, = 
implies self duality and T x i y i = implies anti-self-duality. 

• basis2form(): 

This command finds the basis 2-forms using the definitions in equation 49 of 
|15j . The output can be reached by 
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> LOO; > LOl; > LIO; 

> LtildeOprOpr ; > LtildeOprlpr ; > LtildelprOpr ; 

for L°, Lq, LqI, Lq!, Ly respectively. Here, "tilde" means the variable has 
a tilde and "pr" means "prime" in the meaning of a different spin frame. 

• curv2form(): 

This command calculates the curvature 2-forms using equations 90-91 of [TB] 
and sets them to these variables: 

> ThetaOO; > ThetaOl; > ThetalO; 

> ThetatildeOprOpr ; > ThetatildeOprlpr ; > ThetatildelprOpr ; 

for Oq, 0j, O 1 ] 1 , Qqi, Qqi, Qi> respectively where "tilde" means the variable 
has a tilde and "pr" means "prime" in the meaning of a different spin frame. 

• topologicalnumbers(): 

This command calculates the integrands of the Euler number and the Hirze- 
bruch signature curvature part integrals using the relations 115-117 of [T5] 
namely, 



X 



4tt 2 



/ [|^o I 2 + 4 |^i | 2 + 3^ + |^o I 2 

JM L 

+ 4|^ 1 | 2 + 3^ 2 -2(|$ 00 | 2 + |$ 02 | 2 ) 

- 4 ( I $01 1 2 + I $11 | 2 + I $12 | 2 - 3A 2 ) 1 I A I A m A m (12) 



T 



— f [|^ | 2 + 4|^i| 2 + 3^ 2 - |^o I 2 



Qn 2 jm 

-4|^!| 2 -3^ 2 ] lAlAmArh - 7] s (dM) (13) 

7] s {dM) being the eta-invariant and this value will not be taken into consid- 
eration in the program. The output can be reached by calling 

> eulernumber_integrand; 

> hirzebruch_signature_integrand; 

These numbers have a special importance for the Atiyah-Patodi-Singer index 
theorem of operators on manifolds with boundary |23j [24J [25J [26J . 
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3 Examples 



In this section, we will apply our program to two instanton metrics and 
calculate some objects using the special commands. Lengthy output values 
are not written to avoid distracting the reader's attention. 

3.1 Example 1: Calculations for the Eguchi-Hanson 
metric 

Eguchi-Hanson instanton [T7] is the most similar to the Yang-Mills instanton 
of Belavin et al. [27] and the metric is given as 




Here, 

a x = -(— cos £d9 — sin 9 sin t;d(j)) , 

o y = -(sin £d8 — sin 9 cos £d<fi) , (15) 

a z = -(— d£ — cos 9d<p). 

and the dyad was given in eq. [1] and eq. [2j We run our program in Maple 
and type the name of the input file: eguchihanson 
Now, let us calculate the connection 1-forms: 

npinstanton> connlf orm() ; 

The program calculates and shows the calculated values. Then the gauge is 
found by the program as: 

SELF DUAL GAUGE because all connection 1-forms without 
tilde are zero 

before prompting for another calculation, the procedure shows the output 
variables which can be used for a later calculation, for example, let us call 
of equation 36 of [T5] : 

npinstanton> GAMMA01; 
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To find the Weyl scalars and the Petrov class, one can run, 



npinstanton> weylscalar () ; 

After the Weyl scalars are shown, the Petrov type is found to be: 

Petrov-type D according to anti-self-dual part 

All these results agree with the literature [22J. 
We can find the scalar operator by using 

npinstanton> scalaroperator () 

command. The scalar equation can be solved in terms of hypergeometric 
equations in four dimensional case by the command, 

npinstanton> pdsolve(scalarop, INTEGRATE) ; 

We can easily find the solution for the scalar equation with one extra di- 
mension added trivially if this dimension is a Killing direction (i.e. ds 2 = 
—dt 2 + ds 2 , ds\ being the original instanton metric similar to the case in 
reference [E]). As an example, we just add — k 2 $(r, 6) to the equation after 
dividing by the exponential part which comes from the solution of the Killing 
directions: 

npinstanton> pdsolve(scalarop/ (exp( (m*xi+n*phi) *I) ) 



Here, k t is the eigenvalue corresponding to the extra dimension. We see 
that the addition of the extra dimension results in a more singular equation, 
confluent Heun equation as the solution of the radial part. The occurrence 
of Heun equations in higher dimensional solutions is known in the literature 
(see [19] and references therein) and our example is a new one. 

3.2 Example 2: Calculations for the Nutku helicoid 
metric 

The NP legs of the Nutku helicoid metric can be chosen as [28J, 




-kt*kt*Pni(r,theta) .INTEGRATE) ; 





(16) 
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m M = :(0, 0, cosh(x — z9),i sinh(x — i0)). 

V sinh 2x 



(17) 



This is an example of a multi-center metric. This metric reduces to the flat 
metric if we take a = 0. Since this solution has curvature singularities, it has 
not been studied extensively aside from four articles [18j[19j[20j[29j. 
Let us calculate the basis 2-forms: 

npinstanton> basis2f orm() ; 

After the calculation, to call Lq of equation 49 of [T5] . 

npinstanton> LOO; 

and the output will be 

-\la 2 cosh(x) sinh(a;)(ixl A dx2-^Idx3 A dxA 

dx's were defined previously, at the very beginning of the session in terms of 
the real coordinates. We can calculate the Dirac equation by, 

npinstanton> dirac(); 

and 

npinstanton> dirac [3] ; 

calls the third component of the Dirac equation vector which can be simplified 

and used for calculations. 

We can calculate the connection 1-forms as, 

npinstanton> connlf orm() ; 

For this definition the gauge turns out to be: 

ANTI-SELF DUAL GAUGE because all connection 1-forms with 
tilde are zero 

and the connection 1-forms are shown. 

For the Weyl scalars and the Petrov class, one can run, 

npinstanton> weylscalar () ; 

After the Weyl scalars are shown on the screen, the Petrov type is found to 
be: 
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Petrov-type I according to self -dual part 

The whole results agree with the literature [28]. 

Now let us take a metric of the form ds 2 = —dt 2 + ds 2 , ds 2 being the original 
Nutku helicoid metric. In this case, we need one more 7 matrix for the extra 
dimension (7*). Firstly, let us use our command to find four dimensional 7 
matrices: 

npinstanton> gammamatricesO ; 

It is clear that the extra 7 matrix will be diagonal as the four dimensional 
ones are non-diagonal. Then we can write these commands to form equations 
using {7^, Y} = ■ 

extragamma : =Matrix ( [ [ xl, x2, 0, ] , [ x3, x4, 0, ], 

[0, 0, x5, x6 ],[ 0, 0, x7, x8 ]]): 

We have eight unknowns, so we will have eight equations: 



eqnl : =( 


gamma. 


.contr [1] 


. extras 


ramma+extra^ 


ramma 


gamma. 


.contr [1]) [1,3] ; 


eqn2 : = ( 


gamma. 


.contr [1] 


. extras 


£amma+extra£ 


;amma 


gamma. 


.contr [1]) [1,4] ; 


eqn3 : = ( 


gamma. 


.contr [1] 


. extras 


£amma+extra£ 


;amma 


gamma. 


.contr [1] ) [2,3] ; 


eqn4 : = ( 


gamma. 


.contr [1] 


. extras 


ramma+extragamma 


gamma. 


.contr [1]) [2,4] ; 


eqn5 : = ( 


gamma. 


.contr [2] 


. extras 


£amma+extra£ 


;amma 


gamma. 


.contr [2] ) [1,3] ; 


eqn6 : = ( 


gamma. 


.contr [2] 


. extras 


£amma+extra£ 


;amma 


gamma. 


.contr [2]) [1,4] ; 


eqn7 : = ( 


gamma. 


.contr [2] 


. extras 


£amma+extra£ 


;amma 


gamma. 


.contr [2] ) [2,3] ; 


eqn8 : = ( 


gamma. 


.contr [2] 


. extras 


£amma+extra£ 


;amma 


gamma. 


.contr [2] ) [2,4] ; 



and we solve them, 

sol : =solve (-[eqnl=0 , eqn2=0 , eqn3=0 , eqn4=0 , 
eqn5=0 , eqn6=0 , eqn7=0 , eqn8=0}- 
, {xl , x2 , x3 , x4 , x5 , x6 , x7 , x8}) ; 

Then, we have the form of the unknown matrix: 



preresult : =subs (sol , extragamma) ; 

To satisfy {7*, 7*} = —2, we can set x5 = i,x8 = —%: 

extramatrix : =subs (x5=I , x8=-I , preresult) ; 

Then we can use this matrix with the name "extramatrix". 
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4 Final remarks 



We introduced a Maplell+GRTensorll based symbolic calculator consisting 
of procedures for instanton metrics using a Euclidean Newman-Penrose for- 
malism. A limited version which can run in older versions of Maple is also 
available. 

As the program consists of procedures for different calculations, development 
of the program according to the user's needs is easy by adding procedures. 
The program and sample files can be downloaded from the author's webpage: 

\protect\vrule widthOpt\protect\href {http : / / atlas . cc . itu . edu . tr\string~birkandant/ 
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